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Abstract 


In this paper, a new state of charge estimation method for lithium battery has been presented. Contrary to the conventional methods which 
use complicated battery modeling, a simple resistor-capacitor battery model was used in order to reduce calculation time and system resource. 
Modeling errors caused by the simple model are compensated by the sliding mode observer. The structure of the proposed system is simple, but it 
shows robust control property against modeling errors and uncertainties. The state equation for battery model and the systematic design approach 
for sliding mode observer also have been presented. The convergence of proposed observer has been proved by the Lyapunov inequality equation 
and the performance of system has been verified by the sequence of urban dynamometer driving schedule test. The test results show the proposed 
observer system has robust tracking performance with reduced calculation time under the real driving environments. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Lithium-ion or lithium-polymer batteries are widely used in 
the mobile equipment, electric vehicle, space and aircraft power 
systems for their high energy density, high galvanic potential and 
long lifetimes compared to the lead-acid battery or nickel-metal 
hydride batteries. Due to the complex chemical and physical 
process of the battery, the behavior of battery is hard to pre- 
dict compared with that of electric and mechanic devices. There 
are several parameters associated with battery behavior. The 
key parameter will be the state of charge (SOC) of the bat- 
tery. The SOC corresponds to the stored charge available to do 
work relative to that which is available after the battery has been 
fully charged. SOC can be viewed as a thermodynamic quantity, 
enabling one to assess the potential energy of the system. Since 
the battery charge/discharge current control is based on the SOC 
information, the correct indication of SOC are of considerable 
importance such as hybrid electric vehicle (HEV) application. 
Since there is no sensor available to measure SOC, it should be 
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estimated from physical measurements by some mathematical 
algorithm. 

There have been many attempts to estimate the SOC of bat- 
teries. The most common methods are the charge counting and 
Kalman filter approach. Charge counting or current integra- 
tion is the most commonly used technique, requiring dynamic 
measurement of the battery charge—discharge current, the time 
integral of which is considered to provide a direct indication of 
SOC [1]. However, the charge counting is an open loop SOC 
estimator and thus the errors in the current detector are accu- 
mulated by the estimator. The error is due to noise, resolution, 
and rounding are cumulative and large SOC errors can result. 
A reset or recalibration action is, therefore, required at regu- 
lar intervals—in all electric vehicles. This may be carried out 
during a full charge or conditioning discharge, but is less appro- 
priate for standard HEV operation where full SOC is rarely 
achieved. 

Cell-impedance measurements have also been reported as a 
useful technique for resetting or adjusting SOC estimates from 
integration based methods. However, from results of various 
studies undertaken to identify the impedance variation of battery, 
with SOC, contradictory views to their usefulness in practical 
systems currently remain unresolved [2]. 
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Nomenclature 

A, B, I’, E system parameters 

Ch nominal capacity of the cell (Ah) 

Cp polarization capacitance (F) 

e error states 

Afi, Af, Afs modeling errors and uncertainties 
H gain matrix 

P, Q, Pf, Qf, R positive definite gain matrix 

Rp diffusion resistor (Q) 

R; ohmic resistance (Q) 

Voc(Z) open-circuit voltage as a function of SOC Z 
Vp polarization voltage (V) 

Vp estimate of the polarization voltage (V) 
V; cell terminal voltage (V) 

Ŷ, estimate of the cell terminal voltage (V) 
x system states 

y output state 

Z state of charge 

2 estimate of the state of charge 

Greek letter 

p switching gain 


The Kalman filter methods are well known technology for 
dynamic system state estimation such as target tracking, nav- 
igation and battery field [3,4]. It provides a recursive solu- 
tion to optimal linear filtering, for both state observation and 
prediction problems. The unique advantage of the Kalman 
filter is that it optimally estimates states affected by broad- 
band noise contained within the system bandwidth. The draw- 
backs of Kalman filter are difficulties of feedback gain selec- 
tion. If the gain is not properly selected, the estimated states 
will diverge. Also Kalman filter have some limitations for 
a real implementation such as perfect modeling of the plant 
and Gaussian distribution of the external noise. If these con- 
straints are not satisfied, the performance of the Kalman fil- 
ter will be degraded and thus cannot be used in the real 
applications. 

Other reported methods for estimating the SOC have been 
based on artificial neural networks and fuzzy logic principles 
[5,6]. Since these techniques incur large computation overhead 
on the battery management controller, it can be a problem 
for online implementation. If high performance digital sig- 
nal processing (DSP) chips are used as controller, neural net- 
works, in particular, can reduce the calculation loads for large 
number of empirically derived parameters required by other 
methods. 

In this paper, the new sliding mode observer design method 
has been proposed for the battery SOC estimation. The pro- 
posed sliding mode observer can overcome the above mentioned 
drawbacks by using sliding mode techniques. Sliding mode con- 
trol was first introduced by Utkin [7]. It is a technique robust 
in the presence of parameter uncertainties and disturbances. It 
entails the construction of an equilibrium manifold and a con- 


trol, designed to drive the system state to the sliding manifold 
and maintain it on the equilibrium manifold. The equilibrium 
manifold is constructed so that the system restricted to the man- 
ifold has a desired behavior. As the sliding mode controller has 
been known to have the robustness under the presence of param- 
eter variations and disturbance, the sliding mode observers are 
also known to have similar robustness properties [8—10]. The 
main characteristic of the sliding mode observer is that it has 
robust tracking performance under modeling uncertainties envi- 
ronments and simple control structure. 


2. Battery modeling 


There have been many attempts to develop battery modeling. 
The most common methods are electro-chemistry model and 
electro-circuit model. While detailed chemistry-based models 
have been built to investigate the internal dynamics of the bat- 
tery, these models are generally not suitable for electrical system 
design approach [11]. On the other hand, circuit-based models 
have been built by the electric circuit parameters such as capac- 
itor, resistor, voltage source and so on [12]. It is commonly used 
method for battery controller, because it is possible to express as 
mathematical formulas. It has been known that perfect battery 
modeling is hard to achieve for every operating conditions using 
the electro-circuit model. Therefore, many complicated electri- 
cal modeling methods have been developed to reduce modeling 
errors. However, these methods increase the calculation time, 
system complexity, resources and they can be a cause of insta- 
bility. 

A simple resistor—capacitor model is employed to the lithium 
battery modeling in this paper. All of the modeling errors, uncer- 
tainties, and time varying elements are considered as external 
disturbance at all. The merit of this model is simple, little com- 
putational time, and the modeling errors are compensated in the 
robust sliding mode observer. 

A resistor—capacitor electrical model of lithium-polymer bat- 
tery consists of non-linear voltage source Voc(Z) as a func- 
tion of SOC Z, a capacitance C, to model chemical diffusion 
of the electrolyte within the battery, a diffusion resistance Rp 
as a function of current J, an ohmic resistance R, and termi- 
nal voltage V;. The resistor—capacitor electrical model includ- 
ing uncertainties is shown in Fig. 1. The terminal voltage is 
given as 


Vi = Voc(Z) + IR; + Vp (1) 
The time derivative for SOC Z can be expressed as follows: 


I 1 
=—= Ro,” — Voc(Z) — Vp) (2) 


where Zis the instantaneous current and C, is the nominal capac- 


ity of the cell. The polarization voltage due to the current is 
shown as 


1 I 


-——Vp + = (3) 
RpCp Cp 


Vp = 
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Fig. 1. The resistor—capacitor electrical modeling of lithium battery. 


The time derivative of V; is obtained by assuming d//dt=0, 
and the complete state equation is given as 


l pd 
Ca RO” Q 


V, = Voc(Z) T Vp = 


Ve Voe(Z) + Lpd g j 
RpCp © Ree Cn C RG 

= —41 V; + a1 Vəæ(Z) + byl, 
Ż = a2 V; — a2 Vəæ(Z) — a2 Vp, 


Vp = —a1 Vp + b21, 
y=[1 0 OIV% Z WI (4) 
where aj =I1/(RpCp), a2 =1/R:Ca), bi =(1/Cn) +(1/Cp) 
+(R/Rp_Cp), b2 = 1/Cp. 

This battery model is not linear if Vo, is not linearly pro- 
portional to the SOC Z. Actually Voc is rather piecewise linear 
with respect to Z. The relationship is developed from the cell 
experimental data, where open circuit voltage (OCV) tests are 
performed on successive discharge of the battery, by the appli- 
cation of periodic current discharge. As for the temperature 
variation from +55 to —30°C, the OCV of a lithium-polymer 
battery (LI-PB) varies non-linearly over the battery SOC as can 
be seen in Fig. 2. 

In order to develop piecewise linear model, define 
Voc(Z) =«Z+d for some range of Z and thus « is not a constant 
but varies depending on the Z. d; is defined as Voc linearization 
error caused by the piecewise linearization 
V, = —aV; +aıkZ + bi I + dı = —a,V, +a Z+b1I +d}, 
Z = aV; — ankZ — a2Vp + d2 = aV; — a2 Z — a2 Vp + dd, 


Vp = —a1 Vp + b21 + d3 (5) 


This model is not accurate compared with the real cell data. 
Therefore, the unknown non-linear disturbances terms are added 
to the model to compensate for the modeling errors 


Open Circuit Voltage [V] 


Fig. 2. Open circuit voltage vs. SOC of the lithium-polymer battery. 


V; = -a V; + a1 Z + bil + Afi, 
Z= aV; — a22Z — a2 Vp + Af, 
Vp = —d] Vp +b: + Af (6) 


where Afi, Af2, Af not only represent non-linearities caused by 
linearization error d; and modeling error, but also time-varying 
terms and internal/external disturbances. It can be decomposed 
as 


Afi = I, Af = Iné, Afs = D3 (7) 


where I1, l2, l'3 are known values and & represents unknown 
quantity whose bound is limited. 

Observability of the linear system can be obtained from con- 
struction of the observability of matrix OM 
T 
] 


Om=[C CA CA? (8) 


where C=[1 0 OJ] and A is the system parameter matrix 
given in Eq. (6). When expanded, Om is 

1 0 0 i 
Om = a ail 0 (9) 


at + 41142 aay, + es aan + aja 


Under every operating condition, the observability matrix is 
always full rank. Therefore, the suggested modeling is observ- 
able and thus possible to estimate the internal state of the battery. 


3. Sliding mode observer design 


Consider the following time invariant linear system [13]: 
x= Ax + Bu + TE(x, u) (10) 
y0; d1) 


where x€ R”*! is the state vector, A € R”*”, B € R”*™ is full 
rank, u € R”*! is the control input, C € R”*™ such that CB is 
a non-singular matrix, and ye R”*! is the output, I € R”*! 
and £ € R™*! is the bounded disturbance input, i.e., such that 
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|€|| < A. The controllability matrix (A, B) is completely control- 
lable and (A, C) is completely observable. It is also assumed that 
the input-output system is minimum phase. 

A sliding mode observer for the system Eq. (10) is 


xX = AX + Bu + H(y — Ẹ) + pl sgn(y — J) (12) 
y= CK (13) 
where the gain matrix H and the switching gain p are chosen 


so that the stability of the observer system is preserved. The 
discontinuous feedback input is defined as 


+1, ey>0 
sgn(ey) = (14) 


—-l, ey <0 


The state reconstruction error is defined as e = x — ĉ. Sub- 
tracting Eq. (10) from Eq. (12) gives the dynamical reconstruc- 
tion error system as 


é = Ae — HCe + I&(x, u) — pI sgn(Ce) 
= (A — HC)e + TE(x, u) — pl sgn(Ce) (15) 


éy=y—y=Ce (16) 


The feedforward gain matrix H can be obtained in two ways; 
pole assignment method and LQ method. The LQ method is 
easier to obtain gain matrix H using the Riccati equation as 


AP + PAT — PC'R"'CP=-Q (17) 


where Q, R are arbitrary semi-positive definite and positive def- 
inite matrices, respectively, has a positive definite solution P. 
Then AT — CTH! is stable with 


HT = R!CP (18) 


which is equivalent to the stability of A — HC. In fact H is the 
observer gain matrix for the system Eq. (10). 

By using appropriate Lyapunov equation we select a matrix 
W such that the reconstruction error system is asymptotically 
stable. 

Let Pp be the positive definite solution of the Lyapunov equa- 
tion 


(A — HC) P; + P(A — HC)" = — Qs (19) 
where Qf is an arbitrary positive definite matrix. Set 
TTP = WC (20) 


Then the asymptotically stability of the reconstruction error 
system is guaranteed if W is a positive definite matrix. 
A Lyapunov function candidate for Eq. (15) is 


V(e) = e! Pre (21) 
Then 
Vie) = é! Pre + e! Pre = e'((A — HC) P + P(A — HC)" )e 
+ (ris — pr? sgn(ce)) Pre 
+e™P(rté — pr sgn(ce)) 


= —e" Ore + UI" Preg — pI! Pr esgn(ey)) 

= —e! Ore + 2(WCeé — pWCe sgn(ey)) 

= —e! Ore + 2Wleyé — pey sgn(ey)) 

= —e' Ore + 2Wey(& — psgn(ey)) (22) 

The —eT Qe is always negative and the latter part of the equa- 

tion is 
ey(§ — psgn(ey)) <0 fore, > 0, if p >h; 
ey(& — psgn(ey)) <0 forey <0, ifp>h (23) 


where h is a boundary value for &. 
Then the resultant equation will be 


V(e) <0 (24) 
Therefore 

lim e(t) = 0 (25) 

too 


4. Cell parameter extraction 


The large size Li-PB was used for the test. The cell com- 
prises of a LiMn Oy, cathode, an artificial graphite anode and 
is designed for high power application. It has a nominal capac- 
ity of 5.0 Ah and a nominal voltage of 3.8 V. The dimension 
of the cell is 250mm x 125mm x 5mm and weight of the 
cell is 120g. The thermal chamber and the Nittetsu cycler 
were used as charge—discharge equipments for temperature 
regulations. Nittetsu cycler has 0-5 V voltage measurement 
range and 0-120A current measurement. The cycler’s volt- 
age measurement accuracy is -—-5 mV and its current measure- 
ment accuracy is +200 mA. It also has precision ampere-hour 
counter for direct SOC calculation. True SOC was directly 
obtained from this ampere-hour counter. The test was performed 
with fully charged condition to set the SOC to one. As test 
proceeds, the true SOC was calculated by the ampere-hour 
counter. 

Cell characterization tests were performed to extract cell 
parameters. The parameters are based on the nominal data which 
is obtained from 25°C temperature test result. The modeling 
errors and uncertainties values are obtained from the boundaries 
of the operating temperature which ranges from —30 to 55°C. 
This method can be applied to other kinds of batteries by chang- 
ing the nominal parameters. 

This set of test comprised a sequence of constant current 
discharges for 180s and rests for 3600s. The cell started fully 
charged up to 4.2 V before the test begins. The discharge current 
is 5 A and it corresponds to the 1-C rate of the nominal capac- 
ity. This amounts to 5% decrease of SOC for each period. The 
sampled data is collected every second. The purpose of this test 
is to set the OCV over the entire SOC range and the test result 
was shown in Fig. 2. 

The result comparing the electrical modeling with the cell’s 
test data is shown in Fig. 3. It shows the discharge current, true 
cell voltage, modeling cell voltage and the modeling voltage 
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Fig. 3. Current, voltages for true and model cell, and error. 
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Fig. 4. The polarization resistance Rp. 


error. The shapes of the true cell response and the model out- 
put are similar in general, although many details of the true cell 
response are different. This is mainly due to the non-linear char- 
acteristics of the nominal capacitance Cn of the true cell and also 
to the fact that the values of resistances are changed by the SOC. 
The model parameters are calculated to match with the test data 
and the resultant values are Cn =5 Ah, Cp = 200F, R; =0.001 2, 
Ry = 0.003 Q, and Rp is a non-linear resistance which varies on 
the current. The plot of Rp over current is shown in Fig. 4. 


5. Experimental result for SOC estimation 


This data was used to identify nominal parameters 
for state modeling of the battery. The model parameter 


are given aj =1.667, a2 =0.0589, aj) =2.8339, a22 =0.1001, 
bı =0.00672, b2 = 0.005. The uncertainties terms are determined 
by comparing true cell data with model data to minimize the 
errors. The system parameters are 


—1.667 2.8339 0 
A= | 0.0589 —0.1001 —0.0589 | , 
0 0 —1.667 
0.00672 0.2 
B= o |, af=]ļ|o02]ē 
0.005 1 


Assume é is a bounded random signal satisfying ||&|| <0.1. 
Choosing R = 1, Q = h, then the positive definite solution of Ric- 
cati equation Eq. (17) is 


1.3629 0.9530 —0.0046 
P = | 0.9530 1.0241 —0.0077 
—0.0046 —0.0077 0.2999 
From Eq. (18) H is 
1.3629 
H= | 0.9530 
—0.0046 


Note that the eigenvalues of A — HC are —1.5640 + 0.62233, 
—1.6689. Let Qf = 55. The positive definite matrix of the Lya- 
punov equation Eq. (19) is 


3.0645 2.3942 —0.0179 
Pg= | 2.3942 3.6108 —0.0347 
—0.0179 —0.0347 1.4997 


Let W=1 and choosing p=0.2 as switching gain, then 
the resultant observer system error reduces to zero by 
Eq. (24). 

The configuration of observer system is shown in Fig. 5. The 
cell model parameters are obtained by off-line cell test results 
and the sliding mode observer equations are established by on- 
line for charge/discharge current of the LI-PB as can be seen in 
Fig. 5. The charge/discharge current is applied to the LI-PB and 
sliding mode observer simultaneously. The terminal voltage of 
the LI-PB is measured as output and fed into the sliding mode 
observer to compensate for the errors, and output of the observer 
is the estimated SOC. The controller has been built with infineon 
16-bit microprocessor XC 167-40 MHz. The calculation time for 
one cell including current and voltage measurement is around 
10 ms. Conventional methods for Kalman filters may take 50 ms 
for example. 

The results of SOC estimation using sliding mode observer 
are shown in Fig. 6. The estimated model output is controlled 
with respect to cell terminal voltage with switching ripple, and 
the estimated SOC follow the true SOC although it has some 
deviation at the start/end of rest period. This is caused by the 
discontinuous current and is affected by the abrupt change of 
Rp. In the discontinuous period, the sliding trajectory is away 
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å =A + Bu+H(y-ŷ)+ pl sgn(y- 9) 


Fig. 5. Configuration of the proposed observer system. 


from the sliding surface by the discontinuous function, but 
the trajectory tracks into the sliding surface in a short time. 
The one cycle of Fig. 6 is shown in Fig. 7. The estimated 
output voltage tracks cell voltage with chattering ripples. The 
estimated SOC also tracks true SOC with chattering ripples. 
The average value of the estimated SOC is close to the true 
SOC. This result shows the proposed sliding mode observer 
can track SOC accurately although the cell modeling is not 
accurate. 
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Fig. 6. The results of proposed sliding mode observer. 
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Fig. 7. One cycle plot of Fig. 6. 
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Fig. 8. The results of whole UDDS test. 
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Fig. 9. One UDDS test cycle result for model voltage. 


To verify the performance of the proposed observer at the real 
driving situation, the second test was performed as a sequence 
of 20 urban dynamometer driving schedule (UDDS) cycles. It is 
operated by series of charge—discharge pulses and 5 min rests, 
and spread over the 100-0% SOC range. It can be seen that the 
SOC decreases by about 5% during each UDDS cycle. Fig. 8 
shows the result of overall UDDS cycle current, true cell and 
model cell voltage, their voltage error. The modeling error is less 
than 20 mV for 20-80% of SOC range. For clear view, the one 
cycle of UDDS is shown in Fig. 9. The proposed sliding mode 
observer was applied to the overall UDDS cycle. The resultant 
SOC for whole UDDS cycle has been shown in Fig. 10. The esti- 
mated SOC and error for whole UDDS cycles are shown in the 
figure. The SOC error is bounded to 3% of all the cases. The tra- 
jectory of the estimated SOC and error for the one UDDS cycle 
are shown in Fig. 11 in order to show clear view of the sliding 
mode observer behavior. The trajectories are always confined 
to the true SOC with the chattering value. This chattering can 
be smoothed by saturation function instead of sign function. In 
other way, the average value of the estimated SOC can be close 
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Fig. 10. Estimated SOC and error for UDDS cycles. 
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Fig. 11. One cycle result of the estimated SOC and error. 


the true SOC. In this way, the suggested sliding mode observer 
can be directly applied to the HEV environment with superior 
performance. 


6. Conclusions 


The sliding mode observer design method for LI-PB state 
of charge estimation has been presented in this paper. The sim- 
ple R-C model was used for Li-PB modeling and the modeling 
errors or uncertainties caused by the simple model were com- 
pensated by the proposed sliding mode observer system. The 
systematic design method has been presented and the Lyapunov 
inequality equation has proved convergence of the proposed 
observer. The performance of the proposed system has been 
verified by the UDDS cycle test which is very harsh environ- 
mental test. The SOC error is confined to the acceptable level, 
less than 3% in most cases which is applicable to the real 
environments. 
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